Motional frequency shifts of trapped ions in the Lamb-Dicke regime 



I. LizuaiiJ3 and J. G. Mugilll 

Departamento de Qmmica-Fisica, Universidad del Pais Vasco, Apdo. 644, Bilbao, Spain 

J. Eschnei0 

ICFO, Institut de Ciencies Fotdniques, 08860 Castelldefels, Barcelona, Spain 

First order Doppler effects are usually ignored in laser driven trapped ions when the recoil fre- 
quency is much smaller than the trapping frequency (Lamb-Dicke regime). This means that the 
central, carrier excitation band is supposed to be unaffected by vibronic transitions in which the 
vibrational number changes. While this is strictly true in the Lamb-Dicke limit (infinitely tight 
confinement), the vibronic transitions do play a role in the Lamb-Dicke regime. In this paper we 
quantify the asymptotic behaviour of their effect with respect to the Lamb-Dicke parameter. In 
particular, we give analytical expressions for the frequency shift, "pulling" or "pushing", produced 
in the carrier absorption band by the vibronic transitions both for Rabi and Ramsey schemes. This 
shift is shown to be independent of the initial vibrational state. 

PACS numbers: 03.75.Dg, 06.30.Ft, 39.20.-|-q, 42.50.Vk 



I. INTRODUCTION 



There is currently much interest in laser cooled trapped 
ions because of metrological applications as frequency 
standards, high precision spectroscopy, or the prospects 
of realizing quantum information processing [l| . The ab- 
sorption spectrum of a harmonically trapped (two-level) 
ion consists of a carrier band at the transition frequency 
ujQ and first order Doppler effect generated sidebands, 
equally spaced by the trap frequency lot, see Fig. [T] The 
excitation probability of a given sideband, and thus its 
intensity, depends critically on the so-called Lamb-Dicke 
(LD) parameter 77 = [S-^i/ (2m.wr)]^/^, with fc^ being the 
driving laser wave number. If the LD regime is assumed 
[r] <C 1), the intensity of the fcth red or blue sideband 
scales with 77*^ [J, S, S], k = 1, 2, 3, so the number of 
visible sidebands diminishes by decreasing 77. It is then 
usually argued that in the LD regime the absorption at 
the carrier frequency is free from first order Doppler effect 
B 0) • Of course this is only exact in the strict Lamb- 
Dicke limit, 7; = 0, and for high precision spectroscopy, 
metrology, or quantum information applications, it is 
important to quantify the effect of the sideband tran- 
sitions in the carrier peak, in other words, the asym- 
potic behaviour, as 77 ~ 0, of the frequency shift of the 
carrier peak contaminated by vibronic, also called side- 
band, transitions in which the vibrational state changes.^ 
The inverse effect, in which the sideband is shifted by 
a non-resonant coupling to the carrier, has been previ- 
ously studied in the field of trapped-ion based quantum 
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FIG. I: Excited state probability of a trapped ion after a vr- 
pulse has been applied. The ion is initially in the \g, 2) state. 



computers H, 0|- To get insight and the reference of 
analytical results, we shall examine a simplified one di- 
mensional model neglecting decay from the excited state 
(resolved sideband regime [8,]). The shift dependence on 
the various parameters (duration of the laser pulses, Rabi 
frequency ^r, ujt) will be explicitly obtained making use 
of a dressed state picture and a perturbation theory with 
respect to rj. The cases of Rabi and Ramsey excitations 
will be examined separately since they may be quite dif- 
ferent quantitatively and have different applications as 
we shall see. 
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^ Even though several transitions contribute to a given peak, it 
is named according to the dominant transition: thus we have a 
carrier peak or fc-th sideband peaks. 



A. Notation and Hamiltonian 

We consider a two level ion, with ground (|.g)) and ex- 
cited (|e)) states and transition frequency ujq = ijJe ~ ujg, 
which is harmonically trapped and illuminated by a 



2 



monochromatic laser of frequency ul- In a frame rotating 
with the laser frequency, i.e., in a laser adapted interac- 
tion picture defined by Ha = ^huJLcrz, and in the usual 
(optical) Rotating Wave Approximation (RWA), the ion 
is described by the time independent Hamiltonian [§, 



H = hujT {a^a + 1/2) -^az 



(1) 



where A — uj^ ^ loq is the frequency difference between 
the laser and the internal transition (detuning), — 
|e)(e| - \g){g\, (7+ = \e){g\, cr_ = \g){e\, and a, are 
annihilation and creation operators for the vibrational 
quanta. 

Let us denote by \g,n) (|e,n)) the state of the ion 
in the ground (excited) internal state and in the n*'' 
motional level of the harmonic oscillator. In general 
the Hamiltonian ([T]) will couple internal and motional 
states. The {\g,n), \e,n)} states form the "bare" basis of 
the system, i.e., the eigenstates of the bare hamiltonian 
Hb = H{^B. = 0). The energy levels corresponding to 
the bare states are given by 



HA 



E„ 



(2) 



with En — Tiujt [n + 1/2) being the energies of the har- 
monic oscillator. These bare energy levels are plotted in 
Fig. m (dotted lines) [i[l3] as a function of the detuning. 
They are degenerate when A = ±kujT, k — 0, 1, 2, . . ., 
but the degeneracies are removed and become avoided 
crossings when the laser is turned on, see Fig. [2] (solid 
lines). At these avoided crossings transitions will occur 
between the involved (bare) states, which are nothing 
but the mentioned carrier {k = 0) and sideband (A; > 1) 
transitions [l|. The splitting at each crossing gives the 
coupling strenght of a given transition [lo| . and the dy- 
namics of the system is then governed essentially by the 
reduced 2-dimensional Hamiltonian of the involved levels. 

Apart from these resonant transitions, off-resonant ef- 
fects will also take place since, strictly speaking, the sys- 
tem is not 2-dimensional. In particular, near the atomic 
transition resonance (A ^ 0), there will be a finite proba- 
bility, although small, of exciting higher order sidebands, 
which tends to zero in the LD limit (ry — *■ 0). In this pa- 
per we study how these off-resonant effects behave within 
the LD regime, when r] is made asymtotically small but 
not zero. In particular, we study how these effects affect 
the excited (internal) state probability, shifting the po- 
sition of the central resonance, which is crucial in fields 
such as atom interferometry (Hi] or atomic clocks with 
single trapped ions [l^ . where tiny deviations from the 
Doppler free form of the probability distribution could 
affect the accuracy of the measurements. Possible effects 
for state preparation in quantum information processing 
are also studied. 
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FIG. 2: Bare {Qr = 0, dashed line) and dressed {^Ir/ujt = 
0.3, solid line) energy levels (in arbitrary units) as a function 
of the laser detuning. A not to small LD parameter rj — 0.4 
has been intentionally chosen in order to highlight the higher 
order avoided crossings. 



II. FREQUENCY SHIFT 

In precision spectroscopy experiments, the measured 
quantity is usually the excited (internal) state probability 
Pg, regardless of the vibrational quantum number n. If 
a general state of the trapped ion has the form 



I^W) = i9n{t)\g,n) + e„(t)|e,n)] , 



n=0 



the excited state probability will be given by 



Pe=tr(|^)(V'|e)(e|) = ^Pe,n, 



(3) 



(4) 



Tl = 



where P. 



' kn(^)l is the probability of finding the 
|e,n) state. In principle, the sum is over the infinite 
number of available vibrational quantum states, but it 
can be simplified if the LD regime is assumed. In this 
regime the extension of the ion's wavefunction is much 
smaller than the driving laser wavelength, rj << 1, and 
it is possible to expand the Hamiltonian JT]) in powers of 
V, 

Hld = hujT {a^a+1/2) -^^az 

-I- [{I + iria + irja^)a+ + h.c\ , (5) 

which only couples, in first order, consecutive motional 
states. Then, if the ion is initially in the vibrational level 
riQ, only consecutive levels rip ± 1 will be coupled in a 
first order approximation. In other words, only carrier. 
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first blue and first red sidebands will give appreciable 
contributions to Pe(A ~ 0). It is thus possible to keep 
only the no and ng ± 1 vibrational states and restrict 
our study to the 6-dimensional subspace spanned by the 
{\g, no), |e, no), \g, no ± 1), |e, no ± 1)} bare states. The 
excited state probability (|4]) can then be approximated 
by 



e,no + l 



(6) 



in the LD regime. For all numerical cases examined, we 
have checked that adding further vibrational levels and 
using the Hamiltonian (1) leads to indistingishable re- 
sults with respect to the six-state model if the LD con- 
dition is satisfied. 

For an infinetely narrow trap {rj — > 0), only carrier 
transitions are driven (i.e., transitions in which the vibra- 
tional quantum number is not changed) and the central 
(carrier) peak of the excited state probability is exactly 
at atomic resonance, i.e., at A = 0. The generation of 
blue and red sidebands will affect this distribution shift- 
ing the central maximum by 5, where 5 is the detuning 
that satisfies the maximum condition 



dPe^ 

dA 



A=5 



-{P 



A=5 



(7) 

and defines the "frequency shift" in the following sec- 
tions. This frequency shift can be understood as the er- 
ror in determining the center of the resonance , i.e., the 
position of the maximum excitation. It will be shown 
that the position of this maximum, rather than coincid- 
ing with the line center, varies periodically with the trap 
frequency lot when the sidebands are taken into account. 

In the following sections this shift will be calculated 
in different excitation schemes, such as Rabi excitation 
(a single pulse, which is used in atomic clocks as well 
as quantum logic applications); and Ramsey iterferom- 
etry (two pulses applied in atomic clocks and frequency 
standards). 



III. SINGLE PULSE (RABI) EXCITATION 

If an ion is prepared in \ip{ti)) at an initial time ti, the 
state of the system at a later time t / will be given by 



mf)) 



-iH(tf-ti)/n 



mm 



°-/''(e,n|e„)(e„|5,no) 



(9) 



see an example in Fig. [T] For an infinitely narrow trap 
(?7 = 0), Pe{A) is the well known Rabi pattern (solid line 
of Fig. [T]). For non-zero LD parameters, sidebands are 
generated at integer multiples of the trap frequency lot-, 
(dotted line in Fig. [T|). To obtain analytical expressions 
for these partial probabilities we shall follow the pertur- 
bative approach introduced in [lo|. 



A. Perturbative analysis: "Semidressed" states. 

The perturbative approach in [T^ consists on dividing 
the Hamiltonian in Eq. (O as 



Hld — H; 



SD 



(10) 



with 



hA 



HsD = ?iwT {a^a + 1/2) ^cr^ 



a + a' ) (T4 



(cr+ + cr_) , 
(11) 



where Hsd is a "semi-dressed" Hamiltonian, which de- 
scribes the trapped ion coupled to a laser field, but does 
not account for the coupling between different vibrational 
levels. This coupling is described by the term 1/(77). Note 
that Hld reduces to Hsd in the LD hmit (77 0), and 
V{rj) is a small perturbation of Hsd in the LD regime, 
r? < 1. 

Within this perturbative scheme, dressed states and 
energies of Hld are obtained up to leading order in the 
LD parameter 77 in our 6-dimensional sub-space, see Ap- 
pendix \X[ 



where \ea) (ca) are the a*'' dressed states (energies) of 
the system, i.e., eigenstates (eigenenergies) of H. We will 
consider first the case where a trapped ion is prepared in 
a given state jg, no) at time U = and illuminated by a 
single Rabi laser pulse for a time r. 

The partial probabilities are easily obtained by pro- 
jecting the |e, n) state on the state of the system at time 



B. Excited state probability 

With the expressions of the dressed energies (|A4[) and 
dressed states (|A5I) of Hld, one finds, after some lengthy 
algebra from Eq. ([5]), that the probability of finding the 
ion in the internal excited state after a laser pulse of 
duration r is given, for the three relevant motional levels. 



by 

Pe,no — l 



no- 



, , , ^ fir . (jJtt 

(A — ujt) "COS sm ^ 

/„o . . . fir ujtt 
[i I — Awt j sm — cos — 2" 



nt 



' (2no 



l)sm — 



2 4r22 
sin (wtt — riT/2) sin (tJTT + 



e,no + l 



(rio + 1) 



[uJT + ny 



fir LOj-T 

(A + wt) " cos sin ^ 
{yl + Awt) sin — cos — ^ 



(12) 



where £7 = a/JTI^ + A^ is the effective (detuning depen- 
dent) Rabi frequency. 

These probabilities are different from the ones ob- 
tained if counter rotating terms in Hamiltonian ([S]) are 
neglected after appying a motional or vibrational RWA. 
In this case, instead of a six-dimensinal model, three 2- 
dimensional models are solved to yield 



+k 



fk 
J n 



'sm^4^, (13) 



where ^^.n+fc = r2fl(n|e"'("+"^) 4- k) and 



(14) 



These simplified expressions for the excited state proba- 
bilities give quite different frequency shifts as discussed 
later, and do not add to one exactly at one particular 
value of the detuning. 



C. Rabi frequency shift 

We are interested in the behaviour of near reso- 
nance, i.e., A ~ 0. If only leading terms in A are kept 
and the maximum condition ([7]) is applied to the prob- 
abilities in Eq. (fT^ . it is found that for a weak laser 
("weak" meaning here that a = VIh/ujt ^ 1), the fre- 
quency shift oscillates with the trap frequency lot as 



5{t) k. VtRrfa^ f{VlRT)iiii\u)TT, 



with /(^) being the function 



sin^ 



CsinC 



4 sin^ I 



(15) 



(16) 




FIG. 3: Exact (numerical) (dashed line) frequency shift after 
a Rabi pulse of duration r with a fixed Rabi frequency of 
Q.R — 2Tr X lOOHz and 77 — 0.05. The solid lines represent 
approximate upper and lower bounds for the shift and the 
dotted one the shift obtained if the vibrational RWA is applied 
(using Eq. (|13|l for the probabilities). Usual Paul (RF) traps 
have motional frequencies of a few MHz, but in this plot a 
trap frequency of ujt — 2tt x lOkHz has been considered in 
order to distinguish the fast oscilations. A zoom around the 
TT- pulse is shown in the inset. 



The exact (numerical) frequency shift S is plotted in 
Fig. [3] (dashed line) as a function of the pulse duration 
time. The numerical calculations have been performed 
with the full Hamiltonian ([1]), i. e., to all orders in the 
LD expansion, and with a large basis of bare states (more 
than 6). Here, and in all remaining figures, the numerical 
results or the analytical approximation obtained in the 
LD regime are indistinguishable. 

The upper and lower approximate bounds for the fre- 
quency shift (solid lines in Fig. [3]) are obtained at the 
bounds of the fast oscillating term, i.e., replacing sihlutt 
by ±1 in Eq. p^. 



(17) 



If the applied pulse is a 7r-pulse (r^r = tt/^r), the lead- 
ing order contribution to the shift pS]) vanishes and the 
next order in a has to be considered. Under the 7r-pulse 
condition there is some robustness against the shift er- 
ror, reducing the frequency shift to a pulling effect (i.e., 
a positive shift). 



SiT^)^nR^'a'cos'^.^ 



(18) 



which is not zero (see inset in Fig. [3]) except for the values 
of Qr that make the argument of the cosine a multiple 
of tt/2. 

Remarkably, the general frequency shift is inde- 
pendent of the initial vibrational quantum number ng. 
This follows from the fact that the probability for the 
first red sideband is proportional to the initial motional 
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Ion ur/iTi Tj Qr/2-k (Hz) S/2n (Hz) Reference 



''°Ca+ (729nm) IMHz 0.095 10 - 100 10"^^ - lO"'' [13] 
i^*'Hg+ (282nm) few MHz 0.035 10 - 20 10"^" - 10~^^ [14. 
**Sr+ (674nm) 2.5MHz 0.042 250 - 500 10"° - 10"* [16] 

TABLE I: Rabi one pulse excitation (clocks and frequency standards): for ^^^Hg'^, r) and 5 have beenn calculated with 
ujt/2-k = lOMHz. 



Ion ulit/2h rj Q.r/2ti (kHz) 5/2% (Hz) Reference 



Ba+ (650nm) 50kHz 0.26 1.5 - 15 10"^ - 10^ fl7] 
'*"Ca+ (729nm) 2MHz 0.03 5 10"^ fl8] 

TABLE II: Rabi one pulse excitation: Quantum information and quantum logic. 



state TiQ while the first blue sideband is proportional to 
no + 1, see Eqs. (IT^ . When the maximum condition ([7]) 
is applied, the no's are cancelled. Moreover, the result is 
identical to the shift when the ion is initially in the lowest 
vibrational state. In this case, the frequency shift is just 
due to the first blue sideband (no red sidebands exist) 
but riQ — 0. This particular case can be solved exactly in 
a 4-state model, without a perturbative approach, giving 
the same results, see the Appendix IB] 

Note also that if the vibrational RWA is applied and 
the simplified expressions for the probabilities of the ex- 



cited states (Uni are used to compute the frequency shift, 
quite different results are obtained (dotted line in Fig. 
[3]), with particularly high relative errors near the 7r-pulse 
condition. 

In quantum information applications, the parameters 
a and fiij are usually higher than in frequency standards 
since the speed of the operations is of importance, so that 
the shift of the carrier peak may be much larger. We have 
collected some typical numerical values in Tables I and 
11. 



D. Fidelity for a 7r/2-pulse 

The oscillations of the carrier peak shift with respect 
to U)tt, Eq. (fT5|) . may affect other observables as well. 
As an example we find similar oscillations in the context 
of quantum state preparation. When applying a resonant 
7r/2-pulse to a trapped ion initially in the ground state 
the internal state obtained for rj ~ Q in 

|^,,) = i=(|g)+z|e)). (19) 

The contamination due to the higher order sidebands for 
non-zero rj will make the real internal state differ from 
this ideal state. 

We now define the fidelity T as the probability of de- 
tecting the ideal state (fT9|l . 

T = P,d = tr[|^(T))(V'(T)|V'.d)(^»d|] (20) 

= l^Y.\9n{r)-^en{r)\', (21) 

n=0 



see Eq. (|3]), where the sum is in principle over the infinite 
number of vibrational levels. It is plotted in Fig. |4] as 
a function of a = Qh/lut- The fidelity is unity in the 
"ideal" ry = case but smaller otherwise. This fidelity 
oscillates also with the trap frequency, as it is observed 
in Fig. [31 If a 7r/2-pulse is considered, we may rewrite 
the expression for the shift as 

oi svauTT = sm = sm — . (22) 

a 2a 

The maxima of the sin t?- function are marked with cir- 
cles in the abscissa. 



IV. RAMSEY INTERFEROMETRY 

We may also calculate the frequency shift due to gen- 
eration of higher order sidebands in a Ramsey scheme of 
two separated laser fields [Toj . In these experiments with 
trapped ions, one ion prepared in the \g, ng) state is illu- 
minated with two 7r/2-pulses (t^/2 = 7r/2fifl) separated 
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FIG. 4: Probability of detecting the ideal state \ipid) as a func- 
tion of a = Qu/ujt for different LD parameters. The value 
of Qrt is fixed by the (resonant A = 0) 7r/2-pulse condition. 
The circles shown in the abscissa correspond to the maxima 
of the sin function, i. e., a = . , with n = 0, 1, 2, . . .. 
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FIG. 5: Ramsey interference pattern for a non interaction 
time T = 2t and for different LD parameters. An ion initially 
in the \g, 2) state has been considered. 



by a non-interaction or intermediate time T. The state of 
the system at a time 2r^/2 after the two laser pulses, 
in the same laser-adapted interaction picture used before 
{Hq = \ fiujL(Jz) is given by 



|V;(2r,/2 + T)) = e 

(23) 

where Hb — H{ilfi = 0) is the bare Hamiltonian gov- 
erning the dynamics of the system in the intermediate 
region. A simple generalization of Eq. ^ for two sepa- 
rated laser pulses, gives the probability for the different 



transitions. 



S3 j,k a 



{e,n\ef3){ei3\j,k){j,k\ea){ea\g,no) 



(24) 



with Cj^n being the bare energies corresponding to the 
\j,n) bare states (j = g,e), see Eq. The excited 

state probability distribution will be given again by Eq. 
(|4]) in the general case, which is plotted in Fig. [5l It can 
be shown (Appendix [U)) . that for weak lasers, the central 
maximum is shifted by 



S{T) « nRV^a^ 



cos — - — sm — - — 



ujtT 



(25) 



which is also independent of the initial vibrational quan- 
tum number no and where Tt = 2r^/2 + T is the total 
time of the experiment, see Fig. |6l In the T ^ limit, 
this expression reduces to the one calculated for the Rabi 
case when a 7r-pulse is applied, see Eq. For non- 

zero intermediate times T, the leading order in a in Eq. 
(|25)) may be written as 



d{T) 



■ cos sm ■ 



2 + 2 2 ' 

with approximate upper and lower bounds given by 



5{T) « ± 



2 + UrT ' 



(26) 



(27) 



see again Fig. [6] (solid lines). 



V. DISCUSSION 

We have obtained analytical formulae that quantify the 
motional (sideband) effects in the carrier frequency peak 
of a trapped ion illuminated by a laser in the asymptotic 
Lamb-Dicke regime of tight confinement. Estimates of 
the importance of these effects for current or future ex- 
periments have been provided in Tables |T] and [ill The 
importance of the shift discussed here depends greatly 
on the application and illumination scheme. Three dif- 
ferent situations have been considered: 

(a) In single pulse Rabi interferometry, long laser 
pulses are in principle desired in order to obtain narrow 
transitions, since the transition width is proportional to 
1/r, but this is limited by the stability of the laser and 
by the finite lifetime of the excited state. Typical laser 
pulses are of the order of miliseconds, that is, Rabi fre- 
quencies of tens to hundreds of Hertz if 7r-pulses (max- 
imum excitation) are applied, which gives a frequency 
shift of 10"^ to 10-1"* Hz, see TableU Currently, the most 
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FIG. 6: Exact frequency shift (dashed hne) and approxi- 
mate upper and lower bounds (solid lines) as a function of 
the Rabi frequency after a Ramsey 7r/2-pulse sequence with 
intermediate non-interaction time T — 5r. An ion trapped 
within the LD regime {rj — 0.04) with a motional frequency 
iJT — 2-K X 2MHz has been considered. 



accurate absolute measurement of an optical frequency 
has fractional uncertainty of about 10~^^, but frequency 
standards based on an optical transition in a single stored 
ion have the potential to reach a fractional frequency un- 
certainty approaching 10~^^ [III. This means that the 
frequency shift found here corresponds to fractional er- 
rors of the order of 10^^"' — lO^'^" for typical optical tran- 
sitions, which is far beyond the 10~^* level so that the 
shifts can be neglected in this context in the foreseeable 
future. 

(b) This changes significantly for quantum informa- 
tion applications where fast operations are important and 
therefore the shifts are many orders of magnitude bigger 
even in the Rabi scheme, see Table HIl 

(c) Back to metrology, the shift in the Ramsey scheme 
is more significant than in the Rabi scheme, because the 
illumination times are much shorter and thus the Rabi 
frequencies are correspondingly higher. In recent Ramsey 
experiments with the ^*Sr+ ion at 674nm a trap with 
motional frequency ujt ~ 27rx 2MHz (ry « 0.042) is driven 
by a laser with Rabi frequency ^Ir ^ 2tt x 16kHz, which 
corresponds to laser pulses of several iia Tg']. Different 
intermediate times T are used, ranging from T = r to 
T = lOr. It is clear from Eq. ((27| that the frequency 
shift decreases as the non interaction times T increases. 
With these data, Eq. ([27|) gives frequency shifts of (5 w 
2tt X ImHz for T — t, which corresponds to a fractional 
error of order 10~^*. The effect is therefore small today, 
but relevant for the most accurate experiments in the 
near future. 

Finally, a word is in order concerning the physical na- 
ture and interpretation of the shifts studied here. They 
are obviously associated with motional effects induced by 
the laser on the trapped ion, but they do not reflect en- 



ergy level shifts. Our frequency shifts are defined by the 
carrier peak displacement of the excitation probability. 
This probability is calculated with a linear combination 
of dressed states, as in Eqs. ^ and ([M)) . However, 
note that, while the eigenstates are affected (corrected) 
by the laser coupling of motional states characterized by 
the Lamb-Dicke parameter ?], the energy eigenvalues re- 
main unaffected in first order in 77, see the Appendix A. 
Indeed, the exact calculations of the shift (based on the 
general Hamiltonian (1) and converged with respect to 
the number of levels) are reproduced by the approxima- 
tions in which the eigenenergies remain unchanged, i.e., 
as in zeroth order with respect to rj. The carrier peak 
shifts we have examined may in summary be viewed not 
as the result of energy-level shifts but due to dressed state 
corrections which affect the dynamics anyway. A conse- 
quence is their dependence on the illumination time. 
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APPENDIX A: PERTURBATIVE CORRECTIONS 
TO THE SEMIDRESSED STATES. 

The semidressed Hamiltonian ([TT|) is easily digonal- 
ized, with semidressed (i.e., zeroth order) energies and 
states given by 



,(0) 



,(0) ^ 



nn 



1 



(Al) 

|.g,7i) + |e,n) ) , (A2) 



N± being dimensionless normalization factors given by 



2Q 



2 (n±A) 

R 



(A3) 



The dressed energies and states will be calculated by 
standard (time-independent) perturbation theory, with 
the perturbation given by the coupling term V(j]), see Eq. 
()lip . The matrix elements connectiiig the semidressed 
states are given in the LD regime by [10| 



(yn'Sn,n'-l + Vn' + lSn,n' + lj (1 - ^ss') : 



where s and s' is a shorthand notation representing the 
sign (s,s,'= ±). Perturbation theory provides expres- 
sions for the dressed energies 



,(0) 



(A4) 
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with no linear corrections, since the diagonal terms of the 
matrix elements (|A4I) are zero. The dressed states (up to 
linear terms in 77) are given by 

I \ I (0) \ I 'i-V^B.Vn I (0) , , . 

\f ^] - \/°'> V , iy^RV^ . (0) - 



_^ iririR^n + 1 (0) 



(A6) 



|e"+i,±/ 



±0 ' "'T' 



APPENDIX B: 4-STATE MODEL, EXACT 
SOLUTION 



If the ion is previously cooled down to its ground \g, 0) 
state (e. g., via sideband cooling), the problem becomes 
4-dimensional, since no red sideband will be involved, and 
analytical dressed states can be obtained without using 
the perturbative treatment of Section IIII Al In this 4- 
state model the excited state probability will then read 



(Bl) 



within the LD regime. The expressions for the dressed 
eigenergies are given by 



Tiuji 



(B2) 



where s and s' is a shorthand notation representing a sign 
(s,s' = ±). The (angular) frequencies v± are defined 
by j/± = ^{ujT±^Y +Tf^ji, with 17 = + A2 as 

usual. Near A = 0, these v± are frequencies shifted to 
the blue and red with respect to the trap frequency ujt] 
they correspond to transitions among the dressed levels 
and play an important role in the carrier frequency shift 
as we shall see. The corresponding dressed eigenstates 
can be written as a function of the bare states, 



1 



+ 



s'O- A 

i 

s'Vl - A 
V\e,l)], 



— — (ujT + s'n - svs' ) \g, 0) 
l5,l> 

{lot + s'fl — sj^s') |e, 0) 



(B3) 



with Nss' being normalization factors. (Strictly speak- 
ing, these states are "partially" dressed states in the sense 
that they are eigenstates of a part of the full Hamilto- 
nian.) 

If the ion is assumed initially in the ground \g, 0) state 
and is illuminated by a single laser pulse for a time r. 



the probability of \e,n) is 

^e----/^e,n|e.,.,)(e,,v|<?,0) 



P. 



(B4) 



which may be analitically calculated to give 



Pe.o — 



2n 



COS — COS — — 

V 2 2 



/ LOT + ^ Vj^T lot ^ ^ V^T 

+ sm — sm — — 

\ 2 2 



,(B5) 



-Pel = 



/ y^R 



17 - A . 1/+T + A . V^T 

X I sm 1 sm 



(B6) 



These are "exact" results within the LD and four-level 
approximations. The oscillations in Pgo and P^\ may 
thus be viewed as interferences among the dressed states 
contributions and be characterized by frequencies v±. 

Note also that the expressions (jBSP and (jB6[) are valid 
for lasers of arbitrary intensity. In particular, transi- 
tions to higher order sidebands which in principle are 
off- resonant when A = 0, become important when the 
"Rabi Resonance" condition fi^j = lot is fuUfilled. In 
this case Pei reduces to 



1 

Pel ~ ^ sm' 



(B7) 



which shows that terms which are in principle off- 
resonant lead to resonant effects under certain conditions, 

see also 113, mm [23. 

Expressions (|B5p and (jB6p can be further simplified 
by performing an expansion in power series of the LD 
parameter. To leading order in rj, 



P, 



eO 



Pel 



nt 



V 2n 



A+ sin 



nt 

~2 

{lot + 



77202 



cot ■ 



A_ sin 



{lot - 



1 2 



with A± ~ ^^-^ ■ Pei{t) takes the form of a beating 
oscillation witli a fast frequency lot and a slow frequency 

^R- 

The expressions for the excited state probability sim- 
plify when the duration of the laser pulse is fixed. If a 
TT-pulse is applied (r^r = tt/Qh) we have that 



P.. 



^R 



(B8) 



el 



(^) {A,-A.fcos^--^, (B9) 
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which, near atomic resonance (A ~ 0), can be written as 

PeO ~ 1-^, (BIO) 

With these expressions for the excited state probabiHties, 
the shifted position of the central resonance follows from 
Eq. ([7]): the central maximum in Fig. [T]is pulled to the 
right, to higher frequencies, by 



J(t^) « flfU] a cos 



^ 2 ^TT-^ 



(B12) 



the same result obtained in the general 6-state model 
calculation when a 7r-pulse is applied, see Eq. ([18]). 



APPENDIX C: DERIVATION OF THE 
FREQUENCY SHIFT IN THE RAMSEY CASE 



with TV = no {N = uq + 1) for the red (blue) sideband. 
The presence of the blue and red sidebands will shift the 
position of the central resonance to a position satisfying 
the maximum condition ([7|). This gives a shift of 



6{T) « nRT^' 



COS ■ 



(l-a2)2 



a sm ■ 



2 + 



a cos ■ 



uij'T 



Keeping leading order terms in a if low intensity lasers 
are assumed (a — U,r/ujt <C 1) gives the frequency shift 



S{T) « nRTfa 



a cos 



2 + f7fl,T 
2 ^tTi 



cos sm 



sm 



ujtT 



(C2) 



which is Eq. (^5)) . 



From Eq. (I24p and with the (approximate) dressed en- 
ergies (IA4[) and dressed states (|A5[) obtained in Appendix 
[51 we may calculate the probabilities for the different 
\e,n) states. To leading order in the LD parameter and 
near atomic resonance A ~ 0, they are given by 



p 



±1 



± 



(l-a2)2 
aA 



cos 



(2 + TUr) cos ^ 

LOT V 2 



2 

LUrTt 



a cos ■ 



I ^ T 

r2|j Q,R 



sm ■ 



ujtT\ 
4 



. ujtT\ 
a sm 



a sm ■ 



2 



(CI) 
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